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1 Introduction 



Let us consider the motion of an electron in the single band tight-binding approximation. 
It is well known that in the absence of a dc electric field, the electron moves ballistically 
whereas, when a dc electric field is present, Bloch oscillations prevent a current from be- 
ing set up in the system. It is furthermore expected that if the electron is in contact with a 
thermal environment, the resulting scattering mechanisms will suppress the Bloch oscilla- 
tions and lead to a steady current. This phenomenon can for example be described within 
a semi-classical picture of the motion combined with the relaxation-time approximation 
(see HAM! for example). Alternatively, in open systems theory, one describes the thermal 
environment and its coupling to the electron with a global Hamiltonian, and traces out the 
environment variables to obtain, under suitable additional conditions such as weak cou- 
pling and appropriate scalings, effective dynamical equations for the electron alone. In 
treatments of this type, the environment is often described by a set of oscillators (see for 
example [FZJ) or more generally by a free Bose field. For a review of various approaches 
to transport theory we refer to [B]. There has recently been intensified interest in obtaining 
rigorous results along those lines [DFPl lDFl lCDMI. 

We present here a simple, explicitly solvable, fully quantum mechanical and fully Hamilto- 
nian model of such a particle-environment system within the repeated interaction scheme, 
in which the environment is described by a chain of two-level atoms. We show that a dc 
current is created due to the interaction of the particle with its environment. In addition 
to drifting in the direction of the applied field, the electron diffuses around its mean posi- 
tion. We give a full analysis of the probability distribution for the position of the particle 
in the large time regime (Theorem 13 .11) . We then use a repeated measurement scheme to 
describe the increments of the position and energy observables between time and time t 
(Theorems 13.31 and 13.21) . 

The rest of the paper is organized as follows. In Section [2] we give a detailed description 
of the model under consideration. Our main results are stated precisely in Section [3] We 
describe the effect on the particle of the interaction with a single atom in Section @] and 
the main properties of the repeated interaction dynamics in Section [5] The proofs of the 
main theorems are provided in Sections |6] and |7J 

2 Description of the model 

We consider a spinless particle on the one-dimensional lattice Z and submitted to a con- 
stant external force F > 0. The quantum Hilbert space and Hamiltonian of the particle are 
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H p = r(Z), H P = -A-FX, (2.1) 

where A is the usual discrete nearest neighbor Laplacian and X the lattice position opera- 
tor 

-A = J^(2 \x}(x\ - \x+l){x\ - \x)(x + l\), X = y2x\x)(x\. 

Identifying H p with L 2 (T 1 , d£) via the discrete Fourier transform, we have 

-A = 2(l-cosf), X = i% 

Here T 1 ~ [0, 2%[ is the first Brillouin zone and £ the crystal momentum. Defining the 
translation operator 

T = ^2\x + l)(x\ =e" i€ , 

xez 

we can write H p = 2 — T — T* — FX. We note for later reference that 

[X, T] = T, [X, T*] = —T*, [H p , T] = —FT, [H p , T*} = FT*. (2.2) 

When F — 0, H p has a single band of absolutely continuous spectrum, sp(if p ) = [0, 4], 
and the motion of the particle is described by 

T(t) = e itH pTe~ itHp = T, 

X(t) = e itHp Xe- HHp = X + i(T - T*)t, 

showing its ballistic nature. 

When F ^ 0, we set G = 2F~ X sin£ = iF~ l (T - T*). The commutation relations (1231) 
yield 

H p = e iG (2 - FX)e" iG , (2.3) 

from which it follows that H p has discrete spectrum, sp(if p ) = 2 — FX. This is the well- 
known Wannier-Stark ladder. In the momentum representation, the normalized eigenvec- 
tor ipk to the eigenvalue Ep. — 2 — Fk is given by 

^^) = ^ e i(2^ 1 sin«)-feO. (2 .4) 



In the position representation, we therefore have 



Mx) = f 9 * ^mo-v^m ^ = j k . x (1 



where the J v are Bessel functions. From their asymptotic behavior for large v (see e.g. 
Formula (10.19.1) in nOLBCIH we infer that 

I / „ \ \ k ~ x \ 



ibk(x) ~ — = — — r for \k — x\ — > oo, 

* 1 ' ^2n\k-x\ \F\k-x\) 1 1 

which shows that ipk(x) is sharply localized around x = k. The motion of the particle, 
described by 

T(t) = (j*H pTe -itH p = e -itF T 

itH vvM „ , 4 . (Ft\ . , Ft\ } ( 2 - 5 ) 



X(t) = e^Xe- 1 ^ = X + - sin ^— j sin ^ , 
is now confined by Bloch oscillations. 

In what follows, we let the particle interact with one or many 2-level atoms, each of which 
has a quantum Hilbert space "H a = C 2 which we identify with T_(C), the fermionic Fock 
space over C. The atomic Hamiltonian is given by 

H a = dT(E) = Eb*b, 

where E > is the Bohr frequency of the atom and b* , b are the usual Fermi creation and 
annihilation operators. The later satisfy the canonical anti-commutation relations 

b*b + bb* = 1, b 2 = b* 2 = 0. (2.6) 

The initial state of the two-level atoms will be their equilibrium state at inverse temperature 
(3 described by the density matrix 

p p = Z-V^, Zp = Tr(e~ pH ') = 1 + e"^. (2.7) 

The interaction between the particle and the two-level atom is chosen so that its effect is 
to give a right or left kick to the particle, depending on whether the atom is in its ground 
state or in its excited state. More precisely, we set 

V = ^(|ar+l)(x|®6* + \x}(x+l\ g> b) = Tb* +T*b. (2.8) 

To understand this interaction, note that when F > 0, the translation operator T can be 
thought of as a lowering operator for the particle. Indeed, from (|2.4I) one finds 

Tip k = ipk+i (2.9) 
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Similarly, T* acts as a raising operator. As a result, V describes an exchange of energy 
between the two-level system and the particle. We point out that this interaction is very 
similar to the one which appears in the Jaynes-Cummings Hamiltonian where a two level 
atom interacts with one mode of the electromagnetic field of a cavity (through its electric 
dipole moment and in the rotating wave approximation), see e.g. HCDGl [Dull . Thermal- 
ization of the field through repeated interaction with two-level atoms was proven for the 
Jaynes-Cummings Hamiltonian in [|BP[|. The model treated here is very similar to the one 
studied in |[BP|. except that the spectrum of H p , contrary to the spectrum of the mode of 
the electromagnetic field, is not bounded from below. As a result, the system we treat here 
has no invariant state, as we shall see below. 

The full Hamiltonian of the particle interacting with a single two-level system acts on the 
Hilbert space T-L p <g> H a and is given by 

H = Hp + H a + XV, (2.10) 

where A G K is a coupling constant. As the more explicit formula 

H = 2 - T(l - Xb*) - T*(l - Xb) — FX + Eb*b, 

shows, one can also interpret the coupling to the two-level system as altering the hopping 
matrix elements of the original Hamiltonian. The operator H is easily diagonalized by 
noticing that it commutes with the "number operator" 

N = -e iG Xe~ iG + b*b = + ^, (2.11) 

F E 

which has a two-dimensional eigenspace to each of its eigenvalue n G Z. In particular, 
if E — F then the energy H p + H a is preserved by the full dynamics (which will be 
computed in Section [4]). 

We now turn to the description of the repeated interaction dynamics (see e.g. HAPl lBJM, 
BPJ). We let the particle interact successively, each time during a fixed period r > 0, 
with the elements of a sequence of atoms, i.e., during the time interval [(n — l)r, nr[, the 
particle interacts with the n-th atom and with none of the others. The Hilbert space of the 
atomic reservoir is 

M 
n=l 

where M is the number of atoms and each H a ,n is a copy of H a . The Hilbert space of the 
joint particle+reservoir system is H = H p <S> K en v The full unitary evolution U(t, s) of 
the system is thus described by the Schrodinger equation 

id t U(t,s) = H(t)U(t,s), U(s,s) = I, t,8e[0,Mr], (2.12) 
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with time-dependent Hamiltonian 

M M 

H(t) = H p + Y,H^ n + \J2xn(t)(Tb* n + T*b n ), (2.13) 

n=l n=l 

where Xn is the characteristic function of the interval [(n — l)r, nr[ and H a n , b n , b* n are 
the Hamiltonian, annihilation and creation operators of the n-th atom. We will use the 
following notation 

H n = H p + H a , n + \(Tb* n + T*b n ), H n = H n + ^ if a>fc . (2.14) 

l<fc<M 

Note that H(t) = H n when t e [(n - l)r, nr[. 

We denote by B l {V. p ) the Banach space of trace class operators on the Hilbert space H p . 
Given any density matrix for the particle p p G B 1 (H P ), p p > 0, Trp p = 1, we set the 
initial state of the joint system to 

M 

Po = P P ®p'$ M , PT = ®PP- 

n=l 

After n < M interactions, this state evolves into p nr = U(nr, 0)poU(nr, 0)*. To obtain 
the density matrix p Pi „ r of the particle after these n interactions we take the partial trace 
over the environment, 

p Pt , nT = Tr HmvPnT = Tr ncnv U(nr, 0)(p p ® p® M )U(nr, 0)*. (2.15) 

In fact, for reasons that will become clear later, we shall consider the more general linear 
operator defined on B(H P ) by 

D a , n (A) = Tr Wcnv (/® [pf^ M ) U(nr,0) (a® [p^f M ) U(nr,0)\ 

where a is an arbitrary real parameter. Using the cyclicity of the partial trace w.r.t. atomic 
operators and the fact that pp is invariant under the free atomic dynamics we can replace 
U (nr, 0) by e~ irHn U((n - l)r, 0) in the last formula. It then follows that 

D a , n (A) = Tr Wa (J ® p«)e- i ^( J D a , ri _ 1 (A) ® p 1 ^) e iTH = C^D^A)), 

where 

C a (A) = Tr Wa (J ® p^)e-^(A ® e 1 ^. (2.16) 
6 



We conclude that D a>n (A) = C™{A) and hence 

Tr Hp BC n a (A) = Tr(5 ® [p%]® M )U(nr,0)(A ® [p^] ® M )C/(nr, 0)*. (2.17) 

The operator C a will play a central role in our analysis. For later reference, we describe 
its main properties in the following proposition. 

Proposition 2.1 For any a e R, C a is a completely positive operator on i3 1 (H p ) with 
spectral radius 

cosh f(± - a) BE) 

where p is defined in ( 13.21) . Its adjoint w.r.t. the duality (A\B) = TtAB is the completely 
positive operator on B(H P ) given by 

C>l{B) = Tr Wa (J g) py a )e iTH (B g) (%)e* B . 
Going back to the special case a = 0, we have 

so that the discrete semi-group (£o) ngN acting on the density matrices of H p describes 
the reduced Schrodinger dynamics of the particle. The reduced Heisenberg dynamics is 
obtained by duality: for B E B(H P ) and p p G B 1 ^), 

Tr np BC n ( Pp ) = Tr np C* n (B)p p . 

At this point, the choice of M becomes immaterial and we can consider an arbitrary large 
number of interactions. Given an observable B on H p , we write 

(B) n = TrBC n (p p ), 

for its expectation value at time t = nr. 



3 Results and discussion 

We are now in a position to state our main results on the dynamics of the particle. As will 
be shown in Section HI coupling with a single 2-level atom turns the periodic Bloch oscil- 
lations (12.51) of frequency weioch = F into quasi-periodic motion with the two frequencies 

WRioch and 

u = V /(E-F) 2 + 4A 2 , (3.1) 
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(see Equ. (14.31) below). Repeated interactions with 2-level atoms have a much more dras- 
tic effect. The bounded motion of the particle now becomes diffusive. In terms of the 
parameter 

P=^sin 2 (^)e[0,l], (3.2) 
the motion is characterized by a drift velocity 



v A = v d (E,F) = Ztanh(^), (3.3) 



and a diffusion constant 



D = D(E,F) = |- - ptanh 2 (^ H . (3.4) 
More precisely, the following holds. 

Theorem 3.1 Assume that F > 0, A ^ and uqt £ 2ix7L so that p e]0, 1]. Let the density 
matrix p p G B 1 (7i p ) describe the initial state of the particle and denote by \i n the spectral 
measures of the position observable X in the state p Pi „ T , 

»n(f) = J f(x)d^ n (x) = (f(X)) n . (3.5) 
1 . The Central Limit Theorem ( CLT) holds: For any bounded continuous f on R, 



n— >oo 



lim//1 V5SrJ dMx) = i /We " 



2. 7/Tr (X 2 p P ) < +oo, fa 



1) 

n— >oo 77,7" n— >oo 77,7- 



3. TjfTr (e 7 l x lp p ) < +00 for all 7 > a Large Deviation Principle (LDP) holds in 
the sense that, for any interval J C R, 

lim — log fi n (nJ) = — inf I(x), (3.6) 

n— >oo 77, rceJ 

where I(x) is the Legendre-Fenchel transform of t(r\) = log 9 (— rj/ f3E), i.e., 



I(x) = sup r/x — log 9 
Here the function 9 is defined in (|2.18|> . 



/3£ 
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Note that when E = F, the mobility 

v d /3 sin 2 (At) 

u = hm — = . 

*" f^o F 2r 

and the diffusion constant 



D = /x/T 1 (l - sin 2 (Ar) tanh 2 (^pj ^ , 



satisfy the Einstein relation 



lim D = u/T 1 = ak B T. 

F^O 



The rate function in Part 3 is explicitly given by 

/(*) = 



&E . (R{x)-x\\ , f (1 - p)(R(x) + 1)\ . . , t . 



—x 

+oo otherwise, 



where 

a = : r-TTrm n , R(x) = J x 2 + a 2 (l — re 2 ). 

(1 -p)cosh(/3£/2)' v ; v v v 

It is strictly convex on [—1, 1] and satisfies liy^r) = and I(x) > for x ^ v^r. 

Note that the drift velocity and diffusion constant do not depend on the initial state of 
the particle. The CLT gives us the probability to find the particle at time nr in a region 
of size 0(y/n) around the mean value v^jit, whereas the LDP gives information on this 
probability for a region of size 0(n). To put it differently, it yields information on the 
probability that the particle's mean speed falls asymptotically in an interval of size 0(1). 
Loosely speaking, it says that 

Vn{{n{v d + 5v)t}) ~ e-nHto+Svyr) _ 

The peculiar symmetry e(—j3E — rj) — e(rj) immediately leads to the relation I(x) = 
—fiEx + I(—x) which tells us that 

,. ,. 1 , V>n{n[—v - 8v, —v + 8v]t) 

hm hm — log — — — = —pEv, 

SviO n^oo nr Hn\ n [ v — ov , v + 8v\t) 
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1. e., that negative mean velocities are exponentially less likely than positive ones. The 
reader familiar with recent developments in non-equilibrium statistical mechanics will rec- 
ognize here a kind of fluctuation theorem. Indeed, we shall see that the symmetry of the 
function e(rj) is a direct consequence of time-reversal invariance and, in a sense, a remnant 
of the Evans-Searles (or transient) fluctuation theorem (see (13.71) ). 

We have further studied the statistics of the energy changes of the particle, the environment 
and the whole system. Note that the latter is not expected to vanish, since the Hamiltonian 
is time-dependent, so that total energy is not conserved. 

To study the change in the energy of the atomic reservoir we use the following operational 
procedure. The reservoir being initially in thermal equilibrium at inverse temperature 
f3 and the particle in the state p p , we measure the total energy of the reservoir and the 
particle just before the first interaction and just after the n-th interaction. These successive 
measurements yield the four values E Pt o,E Ptn E sp H p and E enVt o, E env<n 6 spH cnv . It 
will be convenient to express the resulting change in energy in terms of the "entropy like" 
quantities 

ASp^ = f3 (-Ep.n — E p o), AS'envjft = — /3(-Eenv,n — -Eenv,o)) 

where 0* = (3E/F. We denote by F n the joint probability distribution of A5 p , n , AS env , n 
and by E n the corresponding expectation. Note that the quantities E enV! o and E erW)n are 
well defined provided the reservoir contains only a finite number M of atoms and n < M. 
However, under those circumstances, AS enV)n does not depend on M and we can therefore 
consider the reservoir contains an infinite number of atoms. This simple thermodynamical 
limit will always be understood in what follows. 

Remark. When applied to electric charge, or more generally to particle number, the two 
measurement processes described above go under the name full counting statistics (see 
e.g. [ABGK] and references therein). The present application is closer to the approach to 
current fluctuations found in HdRL 

Theorem 3.2 1. P n [A5 , P) „ = AS^ n ] = 1. Hence, in the following we set 

— A C - AC 
n L - lkJ p,n L - 1LJ erw,n- 

2. The cumulant generating function ofAS n is given by 

\ogE n [e aASn ] =n log 0(a), 
where the function 9 is defined in (12.181) . 
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3. Its mean value and variance are 



AS n 



n 



-(3Ev d r, 



(AS n + 0Ev d nrf 



n 



4. The CUT holds: For any bounded continuous function f, 



lim E r ' 

n— >oo 



f AS n + (3Ev d nr \-] f 



f(x)e 



-x 2 /2 



((3E) 2 2Dt. 



dx 



2tt 



5. The sequence (P n ) n£ N satisfies a LDP: For any intervals Jcl, 



lim -logP" 

n^oo n 



AS n 



e J 



n 



inf '00), 



with the rate function <p(s) = sup Qg]R (as — \og9(a)). 
6. It satisfies the transient fluctuation theorem 

pn [ASn = _ s l 

L n J_ gTis 



(3.7) 



Part 1 clearly reflects the fact that the number operator (12.111 ) commutes with H so that 

(3*H P + (3H env = (3E(^ + ^ 
is preserved by the repeated interaction dynamics. 

The particle's drift velocity is v d , and one sees therefore that, as expected, its energy loss 
per unit time equals the work done by F per unit time. Simultaneously, the environment 
gains energy at a rate Ev d : indeed, the particle moves on average v d steps to the right per 
unit time, which corresponds to v d elements of the chain gaining an energy E. This leads 
to an average energy gain or loss of (E — F)v d for the full system. In the special case 
E = F, these rates are equal, and the total system neither looses nor gains energy. This 
is a consequence of the fact (mentioned after Equ. (|2.1 11) ) that the interaction term in the 
Hamiltonian commutes with the free Hamiltonian in this case. In general, the total energy 
is not preserved, which is a reflection of the fact that the Hamiltonian of the total system 
is time-dependent, as is clear from (12.131) . 

Note also that the symmetry 0(1 — a) = 9(a) which leads to the transient fluctuation 
theorem (13.71 ) is evident from Proposition ^. U We shall see in Section 15^21 that it is actually 
a consequence of time-reversal invariance. 



11 



Since 



X, 



we expect a very similar result for the position increment AX n = X n — X obtained from 
a double measurement of X at time t — and t = nr. Indeed, the distribution Q n of AX n 
satisfies the following. 

Theorem 3.3 1. The cumulant generating function of AX n satisfies 
g( V ) = lim -logQ" [e" AX "] = logtf 

n->oc n \ phi 



2. Its mean value and variance are 



lim Q n 

n— K>o 



AXr, 



TIT 



Vd, 



(AX n - v d nrf 



TIT 



2D. 



3. The CUT holds: For any bounded continuous function f, 



lim < 

n— ¥oo 



( AX n - v d nr 
V V2D~ 



>71T 



f(x)e 



-x 2 /2 



dx 



2tt 



4. The sequence 



in&i 



satisfies a LDP: For any intervals Jet, 



lim — k>£ 

n— >oo fl 



AX„ 



e J 



n 



inf I(x). 



— 5v , — v + Sv]\ 



5. It satisfies the asymptotic fluctuation theorem 

(3E(v - 5v) < lim - log L 1 rr 

n^oo n Q n e [ v - Sv, v + 5v]\ 

for v g] — 1, 1[ an J 5f > 0, 5ma/Z enough. 



< (3E(v + 5v), 



Remark. This fluctuation theorem is not of transient (or Evans-Searles) type since it only 
holds in the large time limit. It is not of the stationary (or Gallavotti-Cohen) type either 
since there is no stationary state for the dynamics, as we shall see in Section l53l For more 
details on these fluctuation theorems, we refer the reader to e.g. HJPRURMll . 

Note also the similarity with Theorem 13. 1[ Theorem 13.31 describes the position increment 
of the particle without requiring sharp localization of the particle position at time 0, con- 
trary to Theorem 13. 1[ 
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4 Interaction with a single two level atom 



In this section we investigate the dynamics of the particle interacting with a single atom 
described by the Hamiltonian (12.101) . As already remarked, H can easily be diagonalized 
by exploiting the fact that it commutes with the number operator (|2.1 II) . To get a tractable 
formula for the propagator e ltH , it will be more convenient to consider the unitary operator 

U = (Tb*b + bb*) cos 9 - (Tb* - b) sin 9, 

where 9 is chosen such that 

E - F 2A 
cos(2#) = , sin(2#) = — . 

Using the commutation relations (|2.2I) and (12.61) . one easily shows that 

/ ' /// ' = H p -f l-l'o ( // h - \ J + (4, i ) 



It follows immediately that, for F^O, the spectrum of H is purely discrete, 

E-F uq 
sp(tf) = 2 - FZ + ^— 

Ifujo/F is not an integer, all the eigenvalues are simple. The normalized eigenvectors are 
given by 

|0 fe ,_) = U\i) k ) g> |0) = cos^|^fc) ® |0) - sin0|^fc+i) ® |1), 
|0 fci+ ) = t/|-0&> ® |1) = cos^l^fc+i) ® |1) + sin6»|^ fc ) ® |0), 

where |0) and |1) denote the ground state and the excited state of the atom. Of course the 
situation is completely different when F = 0. The spectrum 

sp(tf) = Q(F - F - a*) + [0,4]) |J Q(F - F + u ) + [0,4]) , 

is then purely absolutely continuous. 

From Equ. (14.11) we get the explicit formula for the propagator, 

e itH = U e it(E-F)/2 e itLa (b*b-l/2) e itH pU *^ ^ 

which reduces the proof of the following result to a straightforward calculation. 
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Lemma 4.1 For any operator A on H p and any t £ R one has 

where A t = e ~ ltH pAe ltHp and 

B,M) = ~ ( ~ 9 sm(cot) - — sin 2 (f) ) (AT* - e~ l3E T*A) , 

By noting that X = 2X eg) p@= the preceding lemma can be used to compute the evolution 
of the position observable 

e itH x e -itH = e itH pXe -itH p 

+ f^ {bb * _ b * b) + 2Xi * E J F \ m + T*6)) sin 2 (^j (4.3) 
-i— (Tb* -T*b) sm(u t). 

We conclude that the coupling to a single atom does not substantially alter the long term 
behavior of the particle. In particular, when F^O, the motion remains bounded. We will 
see in the next section that the situation is very different for repeated interactions with a 
sequence of atoms. 



5 Repeated interaction dynamics 

In this section we study the properties of the operator C a , defined in (12.161) . 

It is clear that a i— > C a is entire analytic as a map from C to the bounded operators on 
13 (Hp). Moreover, since 

Tin p BC a (A) = Tr(B <g) P $)e~ irH (A ® p 1 ^) e irH , 

C a is also bounded as an operator on the Schatten-von Neumann class B P (H P ), for any 
V e [l,oo]. 
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Using the cyclicity of the partial trace w.r.t. operators on H a we can write, for a6l, 
where is a bounded perturbation of H, 

ff(a) _ e -a/3H a /2 He al3H a /2 

= H p + H Si + \(e- a ? E/2 Tb* + e afiE/2 T*b). 
This shows in particular that C a is completely positive for real values of a. 

5.1 Gauge invariance 

Lemma 5.1 The operator C a commutes with the evolution of the non-interacting particle, 
i.e., 

C a (e- M Me"^) = e- itH *C a {A)e ltH *, 
holds for all t,oeR and A G B(H P ). 

Proof. From the fact that H commutes with the number operator (12.1 II) we infer 

H = e itN He~ [tN = e itHp/F e itHa/E He~ itHa/E e~ itHp/F 

so that 

e -itH a /E e -irH e itH a /E = Q itH v / F Q -WH e ~itH p /F _ 

Since H p and H a commute, we also have 

e -UH & /E e -iTH^) e itHe./E _ e itH p /F e -irH^ Q -itHp/F 

and hence 

e itHp C a (e~' ltHp A^ tHp ) e~ itHp 

15 



where we used the fact that H & commutes with A ® pp and the cyclicity of the trace. □ 
Introducing the non-interacting evolution operator 

U{A) =e~ irHp Ae irHp , 

we define the "interaction picture" reduced evolution as C a = C a o U^ 1 . Note that, by 
Lemma I5T1 we have 

C n a = C n a oU n =U n oC n a , 

for any n G N. 

5.2 Time-reversal invariance 

Let us denote by C p the complex conjugation on £ 2 (Z), i.e., (C p ip)(x) = ip(x) and set 

C(A) = c p ac;, 

for A £ B(J-L P ). This anti-linear involution implements time reversal of the particle's 
dynamics. Indeed, since the Hamiltonian H p is real in the position representation, C P H P = 
H p Cp, one has C p e ltHp = e~ ltHp C p for alH e M and in particular 

CoU = U- 1 oC. 

Lemma 5.2 For alia G R one has C* a = C o C\- a o C. 

Remark. As will be clear from its proof, this property of C a is a consequence of the time- 
reversal invariance of the dynamics of the particle coupled to a two-level atom. It implies 
that the spectral radius of C a and L\- a coincide. 

Proof. Setting C a (a |0) + ai|l)) = a^|0) + aT|l) defines a anti-unitary operator on % & 
such that C & H & = HJJ & and hence 

C a e itHa = e-^Ca, C aP/3 = p^C a . 

Since T and T* are real w.r.t. C p and 6 and b* are real w.r.t. C a , the Hamiltonian iJ is real 
w.r.t. C = C p <g> C a and one has Ce itH = e~ itH C. 

From the fact that the partial trace on H a satisfies 

Tt^CAC* = C(Tt h A), 
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for all A e B(H P ® H a ), one deduces that, for A,B e B(U P ), 

C(AC a (B)) = C(Tr n M ® P^)e- irH (S ® pj" )^) 
= Tr n ^C(A ® /^)e- iTH (B ® pJ- Q )e iTH C* 
= Tr Wa (C(A) ® p^ H (C(B) ® p^e"^ 
= Tr Wa (C(5) (S> pJ- a )e-^(C(A) ® p^)e^ 
= C(5)£ 1 _ a oC(A). 

It follows that 



Tr Hp AC a (B) = Tv Hp C(AC a (B)) 



Tr Hp C(B)C 1 - a oC(A) 
Ti np CoC 1 _ a oC(A)B, 



and hence =Co£ 1 _ a oC(A). □ 

5.3 Kraus representation 

Since C a is completely positive for a E R, the same is true of C a . The following result, 
describes the Kraus representation of the latter operator. 

Lemma 5.3 For any a € R and A 6 B(T-L p ), one has 

£ a (A) = e a P E p_T*AT + Po A + e~ al3E p + TAT* , (5.1) 



c 



P- = i + e -^ P; Po = i"P, P+ = l + e -p E P- 



with p defined in Equ. (|3.2I) . 

We note that e a/3E p_ +p + e~ a(3E p + = 9(a), where 0(a) is defined in Proposition ^. 11 so 

C a (I) = £*(/) = 0(a)/. (5.2) 

Consider /C Q = 0(a)~ 1 £ Cf as an operator on B 1 ('H P ). Since JC a is completely positive and 
trace preserving, it has unit spectral radius (see [Sch]). Hence, 0(a) is the spectral radius 
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of C a . This proves Proposition ^. 11 The remark after Lemma [5^21 now explains the origin 
of the symmetry 0(1 — a) = 9(a). 

Let us now set a = in Equ. (15.11) and explore the implications of this expression for the 
dynamics of the particle. If p describes the state of the particle, then T* pT (respectively 
TpT*) represent the same state translated by one lattice spacing to the left (respectively 
right). Note moreover that 

P- + Po + P+ = 1, 

so that the reduced evolution £ = £o °^ consists of a free evolution with the Hamiltonian 
H p , followed by a random translation by ±1 or 0, and with probabilities p± or p . Note 
that the dynamics is trivial if p = 0, i.e., if ujqt = 2nm with m E Z. In that case there is 
no translation and the particle evolves according to H p . This can be seen directly on Equ. 
(14.21 ) by noticing that UH P U* = H p + Fb*b. It follows that the propagator factorizes 

and, up to an inessential phase factor and a renormalization of the atomic Bohr frequency, 
the particle and the two-level system evolve as if they were not coupled. This resembles 
the "Rabi oscillation" phenomenon which appears in the Jaynes-Cummings model for 
matter-radiation interaction. In the following we will avoid this resonance and assume 

Introducing a family of i.i.d. random variables Y n , taking the values ±1 and with proba- 
bility p± and p , and defining S n = Y^j=i ^j'> we can ver y concisely write 

£«(p) = e~ inTHp E [T s "pT~ Sn ] e inrHp . (5.3) 

Accordingly, the study of the dynamics of the system is reduced to that of a classical 
random walk. As a further remark, suppose that the initial state of the particle is invariant 
under the uncoupled dynamics, so that 

P = ^2pk\^k)(^k\. (5.4) 

kez 

Then, using Equ. (HjS), we obtain C%(p) = J2k&Pk n} \i ! k){4 ! k\ with 

(n+l) (n) . (n) (n) /c C x 

Pk = P+Pk-i + PoPk + P-Pk+v ( 5 - 5 ) 

Thus, the set of iJ p -invariant states is invariant under the reduced dynamics and the latter 
reduces to a classical Markov chain on this set. 
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Before turning to the proof of Lemma [5731 let us show that the reduced dynamics has no 
stationary state: there exists no density matrix p on T-L p such that Co(p) = p. Indeed, it 
follows from Lemma I5T1 that the subspaces Jd, d G Z, defined by 

J d ={ p G B\n p ) | e" M Ve itHp = e itd p for all t G R} 
={p G B\H P ) | p = ^p h \il>k){ii>k+d\}, 

are globally invariant under C . Hence, if a state p is stationary, so is its diagonal part 

Po = T,kPk\4>k)(4>k\, where p k = {ip k \pip k ). Equ. (O then writes 

Pk-i - Zpp k + e~ pE p k+1 = 0, 

which implies that p k = a + &e^ Efe for some constants a, b G R. But this contradicts the 
fact that 1 = Trp = J] fc p fc . 

Proof of Lemma [5731 We start with the fact that 



rx a -irH A A-a AtH ZqpZ(l-a)0 a -irH a ~ „ Aril 

p e A® Pa e = p a /3 e A ® P(i- a )p e 

Z /3 



so that, applying Lemma HTTl we get 

p^e-^A ® pj- a e ir " = Za/3 ^ (1 " a)/3 p a/3 (Ai-a)/J,t(A t )6*6 + B ( i- Q )/j, t (A t )6 

+£ (1 _ a)w (A t )*6* + C (1 - aW (A t )bb*) . 
Upon taking the partial trace over "H a , we obtain 

and the result follows from Lemma HTTl with some elementary algebra. □ 



6 Proof of Theorem 3.1 



To complete the proof of our main results, we shall need the following technical lemma. 
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Lemma 6.1 1. If F ^ then, for any 77 G R, one /las 

lim ||e i *- H Pe i ' 7 ' X / v/ *e _i *' ffp — e i7?jSf//v/ *|| = 

t— >oo 

2. If F 7^ then, for any r] G R, ?/zere ex/*?* C,, > 1 smc/? ?/za? 



C- 1 < e ~vX/2^H VeV X e -UH Ve - v X/2 < ^ 



/or a// 1 G 



Proof. 1. By Equ. (12.51) . one has e ltHp Xe ltHp = X + B t , where £? t is a uniformly 
bounded operator valued function of t. Duhamel formula yields 

e itH Pe i V x/Vt e -itH p _ e i V x/Vt = V R 

Vi ' 



where 



Rt = f 1 e ^sMx+B t )/Vi Bte isvX/Vi ds 
Jo 



The claim follows from the fact that R t is also uniformly bounded. 
2. By Equ. (12.31) we can write, for any i],t G R, 

Q(t, V ) = jtH Pe -^X e -itH p jr,X = e iG e -itFX e -iG e -ir,X^G^FX e -iG e if,X_ 

From the commutation relation (12.21) we get, for 9 G R, 

F 

so that 

Q(t,Tf) = e 2i (sinC-sin(C+tF)-sin(5+»7)+sin(5+tF+»7))/F^ 

It follows that 77 !->■ Q(t, 77) G B(% v ) extends to an entire analytic function. Moreover, one 
easily shows that 

C v = sup \\Q{t,r))\\ = e 4sinh (l Im ^)/ F , 



for any 77 G C. 

Now since X t = e ltHp Xe~ ltHp is self-adjoint the subspace V t = Dom(e x<2 ) is dense in "H p 
and such that, for <p G T> t , the vector valued function r] i-> e ir?Xt is entire analytic. For 
•0 G X\, 4> eVq and 7/ G R one has 



( e i^| e i^0) = ^|Q( t)r? )0). 
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n- 2 < p-vX/2 V Xt p -r,X[2 < r 2 



By analytic continuation, both sides of this identity extend to complex values of rj. In 
particular, one has 

(e^ie-^}HVW,^)0}, 

for any r] E R. Since V t is a core of e ±TlXt , the last identity extends to all <p E Dom(e~ vX ) 
and ij) E Dom(e r?Xt ). The modulus of its right hand side being bounded by C- irj \ 
one concludes that Ran(e _,?x ) C Dom(e' ?x *) and 

|| e ^ e -^|| < c iv . 

It follows that 

e - V X/2 eV X te - V X/2 = (e ^/2 e -,X/2 ) * (e ,X t /2 e -^/2 ) < ^ 

and 

which together imply 



□ 



We now turn to the proof of Theorem 13. 1[ First, from Equ. (15.31) we find that 

(f{X)) n = Tl(f(X)m Pp )) = E [Tr [f{T-S^nrH pX ^nrH vT S n) ^] 

Hence, using Equ. (12.51) . we get 

(/(*))„ = E [Tr (/ (T~ Sn (X + B n ) T Sn ) p p )] 

= E[Tr(f(X + S n + B n ) Pp )], (6.1) 

where 

5 n = - sin — sm £ + — . (6.2) 



F 

To prove Part 1, we study the convergence in distribution of the sequence of probability 
measures ji n {J) = p n (V^DnrJ + v^nr), J C R, where //„ is defined in (13.51) . By the 
Levy-Cramer continuity theorem, this is equivalent to the pointwise convergence of the 
sequence of their characteristic functions. We shall therefore prove that for any rj E R 
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Using (16.11) . we can write 



E 
E 



Tr (^j r l{X+S n -v A nT+B n )/^/2DnT p 



:) ir?(5n ~v d nr)/\/2Dm 



The classical CLT implies that the first factor of the last line converges, as n — > oo, to 
e~ v I 2 since, as a simple computation confirms 



— E[S n ] = v d , — Vax[S n ] = 2D, 
nr nr 



(6.4) 



where v d and D are defined in (13.31) and (13.41) . The second factor is controlled by Part 1 of 
Lemma I6T1 which implies 

lim Tr (jv(x + B n )/V2D^ p \ = Um Tr U v x/V2D^ \ 

n— >oo V / n— >oo \ / 

The right hand side of this identity is 1 by the dominated convergence theorem. 
The first assertion in Part 2 follows from 

— (X) n = — Tr(X££(p P )) 
nr nr 

= ^E [Tr((X + S n + B n )p p )] 

= —E[S n ] + —Tr((X + B n ) Pp ), 
nr nr 



since the first term on the last line converges to v d (by Equ. (16.41) ) while the second vanishes 
as n — > oo. The proof of the second assertion is similar. 

To prove Part 3, we shall study the cumulant generating function 

e n (n) = log/i„(e^) = logTr (e" x £"(p p )) . (6.5) 
By Equ. (Ot-dO) and (f23]>, one has 

e n {rf) = logE [Tr (e v(x+Sn+Bn) p p )] = logE [e" 5 "] + logTr (e inTHp e" x e- inT % p ) • 
One easily computes the first term on the right hand side 



logE [e vSn ] = nlogE [e* Yl ] = nlogO 
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Writing the second term as 

and applying Part 2 of Lemma l6.1l we get the bound 

| logTr (e^e^e-^pp) I < ^gC v + logTr (e" x p p ) , 
from which we conclude that 

e{rj) = lim -e n (rj) = log# ( J . 

n->oc n \ phi J 

Since e{rj) is differentiable, the Gartner-Ellis theorem (see e.g. HDZIO implies the LDP 
(13.61) with the rate function I(x) related to e{rj) via the Legendre-Fenchel transform. Fi- 
nally, the symmetry 9(1 — a) = 9(a) translates into 

e{rj) = e(—j3E — rf), 

which implies I(x) = —/3Ex + I(—x). 

7 Full counting statistics 

In this section we start with a precise description of the two measurements processes in- 
volved in the formulation of Theorems I3.2l and l3.3l and proceed then to the proofs of these 
results. 

Suppose that the initial state of the particle is described by the density matrix p p and set 
P = P P ® pf . Let Ai, . . . , A m be commuting self-adjoint operators on the Hilbert space 
H = Hp <E> Henv- We assume the A, to have pure point spectrum. We define the vector 
valued observable A = (Ai, . . . , A m ) and its spectrum sp A = sp A\ x • ■ ■ x sp A m . We 
denote by P a the spectral projection associated to the eigenvalue a E sp A. 

The outcome of a first measurement of the observables A at time t — will be a E sp A 
with probability Tr pP a . After this measurement the state of the combined system is re- 
duced to p' = P a pP a /Ti(pP a ). This state now evolves under the repeated interaction dy- 
namics and, after the n-th interaction, becomes U (nr, Q)p'U(nr, 0)*. A second measure- 
ment of A at time t = nr will yield the result a' with probability Tr U (nr, 0) p'U (nr, 0)*P a i . 
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Thus, the joint probability distribution for the two successive measurements of A is given 
by 

VX(a,a') = TrpPa x Tr U(nr, 0)p'U(nr, 0)*P a / 

= Tr U(nr, 0)P a pP a U(nr, 0)*P a ,. 

Therefore, the probability distribution of the measured increment Aa = a' — a after n 
interactions is 

¥ n A (Aa)= ^U(nr,0)P a pP a U(nr,OyP al . 

a, a 1 £sp A 
a 1 — fl—Aa 

The cumulant generating function of this distribution is 

g n {a) =\og Yl n(A«)e Q ' Aa 

AaSsp A— sp A 

= log Yl ^ a '-^TiU(nT,0)P a pP a U(nr,OyP al , 

a,a'£sp A 

where • denotes the Euclidean scalar product on M. m . At this point, it is useful to note that 

p= Y p «p p «i 

a£sp A 

is a density matrix which commutes with A so that P a pP a = pP a - Hence, we can rewrite 

# n (a)=log Y ^' {a '' a) ^U(nr,0)pP a U(nr,OyP a/ 

a,a'£sp A 

= hgTiU(nT,0)pe~ a - A U(nT,0)*e a - A . (7.1) 
In the special case where the Aj are observables of the particle, we obtain 

g n (a) = \ogTrQ(p p e- a - A )e a - A = logTr p p e'^^e^), (7.2) 
where p p = £ agsp A P a p v P a - 



7.1 Proof of Theorem 3.2 



To prove theorem l3?2l we consider the case A = (l3*H p , —/3H eDV ). From Equ. (17.11) . and 
for any n < M (where M is the number of atoms in the reservoir), we have 

g n (a p ,a cnv ) = logE n [ e apA ^" +aenvASeilv ' ra ] 

= log Tr U(nr, 0)pe- a ^* Hp+acn ^ Hcnv U(nr, o)* e ap/3 ^ p - acnv/3Hcnv , 
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where 

P= Pp,EPpPp,E <8> Pcnv,E'P^ M Peiw^E', 

Eesp Hp 
-B'esp ifenv 

and P Vt E, P C nv,E> are the spectral projections of H p and H env . Since p? M commutes with 
ifenv, this reduces to p = p p <g> pjf M with p p = J2 EespHp P p ,eP p P p ,e- Hence, invoking 
(12.171) . we obtain 

9n{fiipi Ol env) 

= log Tr C/(nr, 0)(p p e " ap/3 * Hp g> p^ M e a ^ Hmv )U(nr, 0)*(e apf} * Hp ® e" —^— ) 

= logTr (e Qp/rHp ® [p^]® M )?7(r2r,0)(p p e- ap/3 * Hp ® [pJ- a °-]® M )C/(nr, 0)* 

= logTr Wp e^^^£" cnv (p p e- ap ^^) 

= log Tr Wp W™(e ap ^ p )£^ cnv (p p e^ B ') 

= logTr w e ap P* Hp/2 £l (e~ ap ^ Hp/2 p p e -<*pP*H P /^ e a P i3*H p /2_ 

Again, the number M of interactions is now irrelevant and we may consider arbitrary 
values of n. The commutation relations (12.21) imply e vHp Te~ vHp = e~ vF T so that, by 
(ED, 

e a p /3*// p /2£^ ( e -a P P*H p /2 A Q -a p f3* H p /2^ ap p* H p /2 

— e (a p +a cnV )l3Ep rp* J^rp _|_ ^ _|_ g - (a p +«cnv )f3Erpj^ rp* 

It follows that 

c/ n (a p ,a en v) = logTr/£^ p+Qcnv (p p ) = logTrp p £^ +Qenv (I) = nlog#(a p + aw), 
from which we conclude that 

logE n [ e «(As P ,„-As cnv ,„)j = n i og ^(o) = o, 

which proves Part 1, from which Part 2, 

g n {a) = logE n [ e aASp <"] = n log 9(a), 

immediately follows. Differentiation of the last identity at a = gives Part 3. Since 9(a) 
is an entire function of a such that 9(0) = 1, log 9(a) is an analytic function of a in a 
complex neighborhood of and Part 4 follows from the Bryc theorem (see fBry|). The 
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Gartner- Ellis theorem directly applies to give Part 5. Finally, Part 6 is a direct consequence 
of the symmetry 0(1 — a) = 9(a) which implies 

log pn i AS P,n = s] e~ s e as = log pn 1 AS p^ = s) e (1 " a)s 

s s 

= g n {\ - a) = g n (a) = log^P" [AS p , n = s] e as , 

s 

and hence P n [A5 Pi „ = -s] e~ s = P™ [AS p>n = s]. 



7.2 Proof of Theorem 3.3 



We now consider the case where A = X. The cumulant generating function of the incre- 
ment AX n is given by Equ. (17.21) . 

g n ( V ) = logQ" [e" AX "] = logTr^ppe-^e"* 

Using the factorization Cq = U n o £™, we further get 

g n {ri) = \ogTrC n (e-^ 2 p p e-^ 2 )U* n ^ x ) 
= logTr(e^/ 2 ££(e-^/ 2 p p e-^^ 

We note that e^ x Te~^ x = e^T and e^ x T*e^ x = e^T* so that, by (TO . 

e^^ofe-^^e-" 1 /^* 2 = e^p^AT + + e v TAT* = Z_ v/ p E {A). 

It follows that 

^(77) = logTr£!! V/3i; (p p )(e-^/ 2 e i ^e'' x e- i "^e-^/ 2 ). 
Part 2 of Lemma I6TT1 yields the estimates 

log C" 1 IV Z% /pE {p p )I < 9n ( V ) < hgC v TrC n _ v/ p E (p p )I, 
and since Tr C n _ v/pE {p p )I = Tr p P C*_ n v/f}E (I) = 6{- V /f3E) n , we finally get 

^) = loge(-^) + o(i), 
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so that 



1 ( V 

g(rj) = lim -g n {rj) = \og9 ' 



rz->oo n \ (3E / 

which proves Part 1 . Part 2 and Part 4 follow from the Gartner- Ellis theorem while Part 3 
follows from the Bryc theorem. Finally, the LDP implies 

.. 1, Q" G [-Q - 8, -q + 5]] . 

hm — log — — - rAy : -y^- = — ml I(— x) + ml I(x) 

«->oo n Q n G [q - 5, q + 5}] \x- q \<8 \x-q\<8 

and the fact that I(—x) = I(x) + (5Ex leads to Part 5. 
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